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SPECIFICATION 

TITLE OF THE INVENTION 
DESIGN METHOD FOR MODE CONTROL DEVICES 

REFERENCE TO RELATED PROVISIONAL PATENT APPLICATION 

This application claims the benefit of United States Provisional Patent 
Application Serial No. 60/010,160 filed January 18, 1996, the entirety of which is 
hereby incorporated herein by reference as if set forth in full herein. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The invention is directed to the field of microwave and optical waveguides. 
More particularly, the invention is directed to methods for designing compact mode 
control and converter elements for use with microwave and optical waveguides. 

2. The Background Art 

Mode control and conversion devices find wide application, for example, in 
microwave heating of plasmas for fusion experiments, multi-mode feeds for RADAR 
systems, microwave waveguide transitions, couplers and mode filters, microwave 
waveguide mode launchers for waveguide transmission systems and for laboratory 
demonstrations, microwave heating systems and optical waveguide couplers, mode 
converters and mode filters optical and microwave gratings, holographic elements, 
filters, phase shifters and the like. 

Mode converters have been extensively used to convert outputs of high 
power microwave sources into lower order modes for plasma heating and low loss 
microwave transmission. 

Currently, periodic gratings are used for conversion of modes in highly 
overmoded circular waveguides. These gratings are formed by periodically varying 
the waveguide radius resulting in a rippled wall structure and are usually analyzed 
by coupled mode theory. Such rippled wall structures result in very slight periodic 
field perturbations. Very high efficiencies have been reported for these gratings, but 
their lengths remain large compared to the waveguide transverse dimension. Various 
techniques have been implemented to reduce the length of these gratings, but the 
overall conversion length remains limited by the grating period. Minimal scattering 
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occurs in such designs and a minimum achievable conversion length appears to be 
equal to one grating period. See, e.g., K. Kumric et al., "Optimization of Mode 
Converters for Generating the Fundamental TE 0I Mode from TE 06 Gyrotron 
Output at 140 GHz/' International Journal of Electronics, Vol. 64 (Jan. 1988), pp. 
5 77 . 94^ and M. J. Buckley, "Compact Quasi- Periodic and Aperiodic TE 0N Mode 
Converters in Overmoded Circular Waveguides for use with Gyrotrons, " I.E.E.E. 
Transactions of Microwave Theory and Techniques, Vol 38 (June, 1990), pp. 712 
- 721 . Similar structures in the form of gratings have been designed for optical 
waveguides. See, e.g., D. Marcuse et aL, "Mode Conversion Caused by Surface 

10 Imperfections of a Dielectric Slab Waveguide, " Bell Systems Technical Journal 
(December, 1969) pp. 3187 - 3215. All of these converter designs are relatively 
lengthy when compared to the radial dimension of the waveguide. A typical 
periodic grating mode converter is diagrammed in FIG. 1. In FIG. 1 a first input 
electromagnetic wavefront 10 is applied to the converter 12. The first wavefront 10 

15 is formed of one or more modes. After interaction with converter 12, which is much 
longer than it is wide, a second electromagnetic wavefront 14 of the selected 
modality is output from converter 12. Such converters are typically 95 - 99.5% 
efficient 

Mode filters for high power microwave sources whose output power is 
20 distributed in various modes, permit extraction of a single mode at the output. 
Previous designs have not proven themselves particularly efficient. See, e.g., /. P. 
Tate, et aL, "Experimental Proof-of- Principal Results on a Mode-Selective Input 
Coupler for Gyrotron Applications", I.E.E.E. Transactions on Microwave Theory 
and Techniques", Vol 42 (Oct. 1994), pp. 1910 - 1917, and U.S. Patent No. 
25 3,771,078 dated November 6, 1973 to H. G. Kidner, et al. 

Mode launchers for efficiently exciting a specific mode into an overmoded 
waveguide can be difficult to construct. However, for a single mode waveguide 
only one mode survives and thus the mode purity is ensured. 

Waveguide adapters, such as tapered sections, are generally used to join 
30 waveguides of unequal radial dimensions. Due to a gradual taper these devices are 
very long as compared to the radial waveguide dimension. See, e.g., W. A. Huting et 
al., "Numerical Solution of the Continuous Waveguide Transition Problem," 
I.E.E.E. Transactions on Microwave Theory and Techniques, Vol. 36 (Nov. 1989), 
pp. 1802 - 1818. 

35 Grating couplers can be used to couple free space light into an optical 

waveguide or vice versa and also for coupling between adjacent waveguides. See, 
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e.g., Nishihara, Haruna and Suhara, "Optical Integrated Circuits", McGraw-Hill 
Optical and Electrooptical Engineering Series, J 989. 

Accordingly, a method for designing more compact, yet equally efficient 
mode converters and control elements for microwave and optical applications would 
be highly desirable. 





WO 98/32187 



PCT/US98/00523 



SUMMARY OF THE INVENTION 



The present invention is directed toward the design of an aperiodic grating 
for mode conversion (defined as any operation on an input set of modes) which is 
based upon an inverse scattering optimization procedure. In accordance with this 
5 method, no predetermined shape of the grating structure is assumed. A constrained 
domain of all surfaces is searched to find an optimum aperiodic conversion surface 
profile for maximum conversion efficiency into the required output mode(s). 
Accordingly, the present method results in the design of rough surfaces or non- 
homogeneous structures for mode conversion. Because this method relies on 

10 scattering produced by short, forceful field perturbations, it is possible to achieve 
very small conversion lengths which are much less than one grating period. 
Precisely because periodic structures are not initially assumed, compact or rough 
structures are generally obtained. This method can be used to design mode 
converters, mode filters, waveguide adapters, mode launchers, power splitters and 

15 combiners, phase shifters and aperiodic gratings for use as couplers in optical 
waveguides. Efficiencies close to the theoretical maximums are obtainable with this 
method. 

The structures achieved by this method have a narrow bandwidth of 
operation and therefore are suitable to the design of filters. 

20 It is known that whenever an obstruction is placed in the path of an 

electromagnetic wave, the energy is scattered into various modes. Such an 
obstruction can be created in a waveguide by varying its dimensions or by changing 
the material of the dielectric inside it. The method of the present invention is 
directed at finding an optimum scattering surface that will scatter essentially all of 

25 the energy in the input mode(s) into one mode or a set of modes desired at the 
output. The following steps are used to design such a surface: 



30 



1. Specify the application of the device (i.e., the frequency of 
operation, type of device, structure, and size of the input and output 
waveguides, the mode composition of the incident electromagnetic 
field and the required mode composition of the output electromagnetic 
field); 



35 



2. Pick a method of variation (i.e., variation of waveguide shape 
only, variation of only the material properties of the obstruction or 
variation of both the shape and the material); 
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3. Choose the material to be used for the obstruction; 

4. Decide the directions in which the obstruction and/or 
5 waveguide shape and material must vary; 

5. Pick a suitable basis function to represent the waveguide shape 
and/or material properties in each of the waveguide dimensions; 

10 6. Choose an initial structure approximation (i.e., give some 

arbitrary (but realistic) values to all the variables (coefficients) in the 
series representations and pick a length, L, for the structure); 



7. Formulate the forward solution; 

15 

8. Perform a global optimization to find values of the coefficients 
that maximize the output power in the desired output mode(s); 

9. (Optional) Obtain multiple designs (i.e., repeat steps 5 through 8 
20 to come up with a number of designs and select an optimal design 

under the circumstances); and 



10. Fabricate the device. 

25 OBJECTS AND ADVANTAGES OF THE INVENTION 

Accordingly, it is an object and advantage of the present invention to 
provide a method for designing compact, efficient waveguide mode control and 
conversion devices useable in microwave and optical applications. 

It is a further object and advantage of the present invention to provide 
30 a method for designing mode converters, mode filters, waveguide adapters, mode 
launchers, power splitters, phase shifters and aperiodic gratings for use as couplers in 
waveguides. 

It is a further object and advantage of the present invention to provide 
a methodology for creating aperiodic mode control and mode conversion surfaces 
35 which are more compact than conventional periodic structures used for these 
purposes. 
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Yet another object and advantage of the present invention is to 
provide a method for the design of narrowband filters, specifically for optical 
applications. 

These and many other objects and advantages of the present invention 
5 will become apparent to those of ordinary skill in the art from a consideration of the 
drawings and ensuing description of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a schematic diagram of a conventional periodic grating used in mode 
10 control and conversion in accordance with the prior art. 

FIG. 2 is a schematic diagram of an input waveguide section, a mode 
conversion/control element and an output waveguide section in accordance with a 
presently preferred embodiment of the present invention. 

15 

FIG. 3 is a schematic diagram of an input waveguide section, a mode 
conversion/control element and an output waveguide section in accordance with a 
presently preferred embodiment of the present invention. 

20 FIG. 4 is a schematic diagram of a mode conversion/control element disposed 

within a waveguide section in accordance with a presently preferred embodiment of 
the present invention. 

FIG. 5 is a schematic diagram of a mode conversion/control element formed of 
25 slabs of materials of varying indices of refraction disposed in an optical structure in 
accordance with a presently preferred embodiment of the present invention. 

FIG. 5A is a schematic diagram of a mode conversion/control element formed 
of slabs of materials of varying indices of refraction. 

30 

FIG. 6 is a diagram showing the construction of a mode converter in 
accordance with a presently preferred embodiment of the present invention. 

FIG. 7 is an exploded diagram showing the construction of a mode converter 
35 in accordance with a presently preferred embodiment of the present invention. 
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FIG. 8 is a diagram showing an assembled version of the mode converter in 
accordance with FIG. 7. 

FIGS. 9 and 10 depict an SMA wave launcher in accordance with a presently 
5 preferred embodiment of the present invention. 

FIG. 11 is a diagram illustrating a simplified scattering structure within a 
parallel plate waveguide in accordance with a presently preferred embodiment of the 
present invention. 



FIG. 12 is a diagram illustrating an optimization surface for TEi to TE2 
conversion in accordance with the scatterer of FIG. 11. 

FIG. 13 is a diagram illustrating an optimization surface for TE t to TE 6 
15 conversion in accordance with the scatterer of FIG. 11. 

FIG. 14 is a diagram showing the construction of a TE! to TE2 converter in a 
parallel plate waveguide in accordance with a presently preferred embodiment of the 
present invention. 



FIG. 15 is a series of four diagrams showing the sequence of optimization of 
the structure for a TEi to TE 2 mode converter structure in a parallel plate waveguide 

in accordance with a presently preferred embodiment of the present invention. 



converter structure in a parallel plate waveguide in accordance with a presently 
preferred embodiment of the present invention. 

FIG. 17 is a chart comparing performance parameters of the three distinct 
30 designs shown in FIG. 16. 

FIG. 18 is a diagram showing the structure of a TEj to TE2 mode converter 

structure in a parallel plate waveguide in accordance with a presently preferred 
embodiment of the present invention. 



10 



20 



25 



FIG. 16 is a diagram comparing three distinct designs for a TE^ to TE 2 mode 



35 
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FIG. 19 is a diagram showing the structure of a TE 02 to TE 0 i mode converter 
structure for circular waveguide in accordance with a presently preferred 
embodiment of the present invention. 

5 FIG. 20 is a plot of efficiency versus frequency for the TE 02 to TE 0 i mode 

converter structure of FIG. 19. 

FIG. 21 A is a plot of width profile versus distance along the z direction for a 
prior art circular waveguide TE02 to TEqi mode converter in accordance with M. 7. 
10 Buckley, et al., "A Single Period TE 02 to TE 0I mode converter in a highly 

overmoded circular waveguide," LE.E.E. Trans. MTT, Vol 39, No. 8 (Aug. 1991), 
pp. 1301 - 1306. 

FIG. 21B is a plot of efficiency versus frequency (GHz) for a prior art circular 
15 waveguide TE 0 2 to TE 01 mode converter in accordance with M. J. Buckley, et ai, "A 
Single Period TE 02 to TEqj mode converter in a highly overmoded circular 
waveguide, " LE.E.E. Trans. MTT, Vol 39, No. 8 (Aug. 1991), pp. 1301 - 1306. 

FIG. 22 is a chart comparing the length and conversion efficiency of the 
20 "Previous" design of FIGS. 21 A and 21B with the "Present" design of FIGS. 19 and 
20. 

FIG. 23 is a diagram of a grating design for a circular waveguide TE^ to TE 0 i 
mode converter for operation at 140 GHz in accordance with a presently preferred 
25 embodiment of the present invention. 

FIG. 24 is a plot of efficiency versus frequency (GHz) for a prior art circular 
waveguide TEq6 to TEqi mode converter for operation at 140 GHz in accordance 

with K. Kumric et al., "Optimization of Mode Converters for Generating the 
30 Fundamental TEqi Mode from TE 0 s Gyrotron Output at 140 GHz," International 

Journal of Electronics, Vol 64 (Jan. 1988), pp. 77 - 94. 

FIG. 25 is a chart comparing the length and conversion efficiency of the 
"Previous" design of FIG. 24 with the "Present" design of FIG. 23. 



35 
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FIG. 26 is a pair of diagrams comparing a pair of designs for a circular 
waveguide TEn to TMn mode converter structure each designed in accordance 

with a presently preferred embodiment of the present invention. 

5 FIG. 27 is a chart comparing the length and conversion efficiency and other 

performance parameters of the "Design 1" design of FIG. 26 with the "Design 2 n 
design of FIG. 26. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

10 Those of ordinary skill in the art will realize that the following description of 

the present invention is illustrative only and is not intended to be in any way 
limiting. Other embodiments of the invention will readily suggest themselves to such 
skilled persons from an examination of the within disclosure. 

The present invention is directed to a new method for the design and 

15 fabrication of mode control devices which can result in aperiodic, totally rough 
structures. These devices can be developed for free space wave propagation or for 
microwave and optical waveguides. Previously it was believed that there are only 
certain type of shapes (mostly periodic structures) that can be used to develop mode 
control devices. In accordance with the method set forth below, however, the 

20 domain of all possible shapes and materials is explored for a scattering obstruction, 
which when placed in the path of the incoming electromagnetic wave will convert 
power from a set of modes in the incoming electromagnetic wave (one or more) to 
any of a set of modes at the output. The method taught here is completely general 
and does not restrict itself to a specific type of pattern variation. Structures obtained 

25 with this method have aperiodic variations and can be very irregular in shape. These 
structures are generally very short as compared to their periodic counter-parts 
obtained using coupled mode theory. 

OVERVIEW 

30 The basic concept is to place an obstruction in the path of the incoming 

electromagnetic field. The shape, size, and material of the obstruction is totally 
arbitrary. When the incoming field interacts with the obstruction, it scatters into 
various modes of propagation. A procedure is then adopted to optimize the shape, 
size, and material of the obstruction so as to maximize the power scattered into the 

35 required mode at the output, so that close to 100% of the input power is transformed 
into the desired mode or that the output mode(s) has the desired characteristic. The 

-9- 
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size and shape of the input and output waveguides is arbitrary. Also the obstruction 
can protrude outside the dimension of the waveguide and the waveguide wall 
dimensions can vary in the mode converter section. 

5 METHODOLOGY 

The variation in the shape of the waveguide and/or material of the 
obstruction, in each of the three dimensions (the three dimensions for rectangular 
and parallel plate waveguides are the Cartesian coordinates: x» y and z; for circular 
waveguides, these are the cylindrical coordinates: r, n and z), can be represented 
10 by suitably weighting a set of basis functions, e.g., delta functions, step functions, 
trigonometric functions, Bessel functions and other functions as would be readily 
understood by those of ordinary skill in the art. 

To illustrate the idea consider an obstruction placed in a circular 
15 waveguide. In the most general case, both the shape of the waveguide and material 
properties of the obstruction will be changed in all directions as shown in FIG. 3. 
However, to keep it simple assume that the shape and material properties of the 
obstruction do not vary in the radial, r, and circular, n , dimensions as in FIG. 2. This 
means that the obstruction 16 fills the whole waveguide 18 and its dimension and 
20 material properties vary only as a function of z (the distance along the waveguide in 
the direction travelled by the incoming electromagnetic wavefront). Such a 
structure would be able to convert only the radial suffix of the incoming mode. 
Specifically, the following conversions are possible: 

TE mn toTE mp 

25 TE^toTM^ 

TM^ to TEjnp 

TMmn to TM mp 

Before the design procedure starts it is necessary to specify the composition 
of the input electromagnetic field and the output electromagnetic field. It is also 
30 necessary to specify some parameters of the output magnetic field which are to be 
optimized to obtain the required output. In the present case, the input is a 
combination of TE mn modes while output is a TE mp mode. The optimization 

parameter is thus the power in mode TE^. 
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The z- variations in the shape and electrical/magnetic material properties of the 
obstruction can be represented by using some suitable basis functions. Some 
examples are given below: 

USING THE WAVEGUIDE SHAPE 
Ar(z) = r 0 n (z) + ^n (z-Z]) + r 2 n (z-z 2 ) + ... (n is the Kronecker delta function) 
OR 

Ar(z) = r 0 U(z) + rjUCz-z^ + r 2 U(z-z 2 ) + ... (U is the step function) 
OR 

Ar(z) = rosin(z) + r!sin(3z) + r2sin(5z) + ... 
OR 

Ar(z) = tqJ 0 (z) + r^z) + r 2 J 2 (z) + ... (J N are the Bessel functions) 
where 0 < z < L ; L is the maximum length of the obstruction. 

If r is the radius of the input waveguide then Ar(z) is the delta change in 
waveguide radius as a function of z. 
Lower Limit on Ar(z): Ar(z) > -r 

Upper Limit on Ar(z): Practical constraint on the largest radial dimension for the 

waveguide. 

USING THE ELECTRICAL PROPERTIES OF THE 
OBSTRUCTION (DIELECTRIC CONSTANT, n ) 

An (z) = e 0 n (z) + e^ (z-ti) + e 2 n (z-t^H-... (n is the Kronecker delta 

function) 

OR 

An (z) = e 0 U(z) + e,U(z-t!) + e 2 U(z-t 2 ) + ... (U is the step function) 

OR 

An (z) = eoSin(z) + e!Sin(3z) + e 2 Sin(5z) + ... 

OR 

An (z) = eoJ 0 (z) + e 1 J 1 (z) + e 2 J 2 (z) + ... (J N are the Bessel functions) 
where 0 < z < L ; L is the maximum length of the obstruction. 
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If n is the dielectric constant of the material in the input waveguide then 
An (z) is the delta change in n as a function of z. 

USING THE MAGNETIC PROPERTIES OF THE 
5 OBSTRUCTION (MAGNETIC PERMEABILITY, ^i) 

A|l(z) = uon (z) + ujn (z-d|) + u 2 n (z-d 2 ) + ... (n is the Kronecker 

delta function) 

OR 

10 A|i(z) = u 0 U(z) + uxUCz-d]) + u 2 U(z-d 2 ) + ... (U is the step function) 

OR 

A|ll(z) = u 0 sin(z) + U!sin(3z) + u 2 sin(5z) + ... 

OR 

A|i(z) = uqJqCz) + uiJi(z) + u 2 J 2 (z) + ... (J N are the Bessel functions) 

15 

where 0 < z < L ; L is the maximum length of the obstruction. 

If \i is the magnetic permeability of the material in the input waveguide then 
A|i(z) is the delta change in fi as a function of z. 
20 Note that for delta functions if the difference between two consecutive Zj/s is 

negligibly small and the difference between two consecutive Ar(z) values is much 
less than a wavelength, then the above representation can be used to model 
continuous surfaces, which is the most general case. For step functions the same can 
be done by restricting the difference between two consecutive Ax(z) and two 

25 consecutive z n 's to be much less than a wavelength. The other two functions 
mentioned above automatically result in continuous surfaces. 

If one type of function is used for the waveguide shape, it is not necessary 
that the same type of function be used for material properties. The functional 
representation of Ar(z), An (z), and A|X(z) can be totally different. 

30 The variables r 0 , r lt r 2 ,...; eo, e lt e 2 , e 3 ,...; uq, u lf u 2 , u 3 ,.„; z h z 2 , z 3 t lt t 2 , t 3 „. M 

d lt d 2 , d 3 ,... for step, delta, sine, Bessel or other functions; and length L uniquely 
represent the shape and material of the obstruction. The power scattered by this 
structure in a particular mode, P 0 , at the output of the structure is then a function of 
all these variables, i.e., 

35 

Pq = Po( r O» r l» r 2""i ^ e l* ^ e 3— "» "0* u l* u 2' u 3' "» z l' ^ z 3>—* etc » L) 



"(2- 
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In accordance with the design methodology of the present invention, one 
optimizes the power P 0 as a function of these variables so that nearly 100% of the 
input power is, for example, converted to the required mode. Optimized values of 
5 the. variables then represent the required design of the mode control device. 

It may be noted here that there is a potentially infinite number of variables 
involved in representing the obstruction in the most general case. To make the 
solution practically viable it is preferable to truncate each infinite series of basis 
functions to a finite number of terms. Thus, one can express Pq as: 

10 

P 0 ( r O r b r 2— > r n> e l> ^ ^ ^m* "0* u l> u 2> u 3 "q* L ) 

where the series for Ar(z) is truncated to n+1 terms, for An (z) it is m+1 terms, and 
A|i(z) is q+1 terms. 

15 Even after truncation, the number of variables for optimization may be very 

large. To optimize say 50 variables simultaneously is often a very difficult problem. 
Most global optimization techniques work best when the number of variables is 
limited to between two and five variables. In order to optimize the mode converter a 
hierarchical, sectional optimization procedure has been developed. 

20 

HIERARCHICAL OPTIMIZATION 
In this procedure one starts with a very course estimate for the features of the 
mode converter. This is done by just a few terms in the basis function representation 
of shape/material and a small length L. This way the initial number of variables is 

25 small. These variables are then optimized in order to obtain a coarse estimate of the 
final structure. Next, more terms are added to the series and L is also increased, thus 
increasing the number of variables and also refining the representation of the mode 
transducing obstruction. Optimization is again carried out and the procedure is 
repeated until a suitable mode converter is identified. Because of a small number of 

30 variables in earlier stages a lot of computation time is saved. 

SECTIONAL OPTIMIZATION 
As one moves up in the hierarchical optimization procedure the task of 
optimization becomes even more difficult due to the increase in the number of 
35 variables. In order to perform a global optimization of all of these variables a 
sectional optimization procedure is preferably employed. For example, assume that 
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one is using a global optimization technique which converges best when the 
number of variables is restricted to two, but the number of variables to be optimized 
is 50. The variables may then be optimized two at a time, starting with the first two. 
While optimizing a set of two the rest of the variables are kept constant. When all 
5 the 50 variables have been optimized once, go back to the first set and repeat the 
sequence. The sectional optimization sequence is continued until a global solution 
is achieved for all the 50 variables. 

It may be observed that as one truncates the series to a small number of terms, 
at the start of hierarchical optimization, one restricts the domain of the structures that 

10 can be formulated for mode conversion. During optimization one searches in this 
restricted domain of structures. If a desired solution is obtained, the process is halted 
and the solution is accepted, otherwise the number of terms in the series is increased 
and/or the length, L, is increased and a search is carried out in a larger domain of 
possible structures. This procedure is continued until an appropriate structure is 

15 identified. Thus no presumption is made about the final shape of the mode 
converter. 

The hierarchical procedure defined above does not result in unique solutions 
for mode conversion structures because the number of free variables is much greater 
than the cost function constraints. Depending upon the initial guess for 

20 optimization, type of basis functions, and number of terms in the series various 
solutions can be obtained for the same problem. Out of the various solutions found 
through the design process only those structures may be considered that are feasible 
to fabricate using the available fabrication procedures. Normally a number of 
solutions would be obtained. Each would be assessed for its feasibility of 

25 construction and the easiest to fabricate would be chosen for development 



30 1. SPECIFY THE APPLICATION OF THE DEVICE 

At this point one must determine and specify the basic known parameters of 
the device to be designed. This requires the following to be determined: 

(a) Frequency of operation, type, structure, and size of the input and 
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output waveguides. (Microwave and lower/upper frequencies: free 
space, parallel plate, rectangular, circular, elliptical, etc. waveguides, 
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and transmission lines like microstrip, stripline, coaxial etc. Optical and 
lower/upper frequencies: free space, slab, fiber etc. waveguides); 

(b) Mode composition of the incident electromagnetic field; 

(c) Required mode composition and characteristics of the output 
electromagnetic field; and 

(d) Parameters of the output electromagnetic field which uniquely specify 
the field and are to be optimized. 

2. PICK A METHOD OF VARIATION 

Pick any of the following methods for generating a wave obstructing 
structure: 

(a) Variation of waveguide shape only; 

(b) Variation of only the material properties of the obstruction; and 

(c) Variation of both the shape and the material. 

Which of the above three methods one chooses depends upon the suitability 
of fabrication for a particular application and waveguide type. Methods 2a, 2b and 
2c may be used for both microwave waveguides and optical waveguides. For free' 
space application there can be two cases: reflection and transmission. In the 
Reflection Case the obstruction would be a rough aperiodic reflecting surface; in 
the Transmission Case a dielectric obstruction with varying shape and properties can 
be used. 

3 . CHOOSE THE MATERIAL 

Choose the material to be used for the obstruction. If method 2a is used then 
the material would preferably be the same as that used for the input waveguide. 

4 . DECIDE THE DIRECTIONS OF PERTURBATION 

Depending upon the input mode structure and the required output mode, 
decide the directions in which the obstruction and/or waveguide shape and material 
must vary. For example for a TE mn mode, in a rectangular waveguide, the following 

conversions can be achieved: 

by varying the structure in y dimension, 
by varying the structure in x dimension. 
15- 



TEmn tO TEmp 
TEmn to TEp n 
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where the "m" subscript refers to the x-direction and the "n" subscript refers to the y- 
direction. 

5 5 . PICK A SUITABLE BASIS FUNCTION 

Pick a basis function to represent the waveguide shape and/or material 
properties in each of the waveguide dimensions. Again, the selection of the basis 
functions would depend upon the particular application and fabrication procedure, 
e.g., sinusoidal basis functions may be more suitable when the mode converter is 
10 going to be pre-formed or molded. By selecting step functions for a circular 
microwave waveguide, a very easy fabrication technique may be implemented as 
shown below. Truncate the infinite series representation to a finite number of terms. 

6. CHOOSE AN INITIAL STRUCTURE APPROXIMATION 

15 Give some arbitrary (but realistic) values to all the variables (coefficients) in 

the series representations. Also pick a length, L, for the structure. Initially choosing 
these values defines an initial "guess" for the shape and material of the obstruction. 
This guess will act as the initial point for the optimization procedure. After some 
experience in dealing with a specific application a designer can use his/her intuitive 

20 feeling as to what type of initial guess would result in faster convergence of 
optimization and best results. See the discussion of initial structure approximation, 
below. 

7. FORMULATE THE FORWARD SOLUTION 

25 Given the information in steps 1 through 5, the forward solution to the 

scattering problem can be formulated. The goal is to find the electromagnetic field at 
the exit plane of the obstruction (O-O 1 in FIG. 1) when a given field is incident on 
the input plane (I-F in FIG. 1). Any boundary value solution technique can be used 
to solve for the output field. A very useful commercial software package for this 

30 purpose is HFSS (High Frequency Structure Simulator) available from the Hewlett- 
Packard Corporation of Cupertino, California. Other software packages are also 
available to perform functions like this as is well known to those of ordinary skill in 
the art. From the solution of the electromagnetic field, the power coupled into the 
required output mode(s), Po, in the output waveguide can be evaluated. Pq is a 

35 function of all the defining variables for the obstruction and needs to be optimized 
(maximized) to find the required structure. 
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8. PERFORM A GLOBAL OPTIMIZATION 

Select a global optimization technique for constrained optimization. It is 
preferable to use constrained optimization because most of the time there would be 
5 realistic constraints on the minimum and maximum values of various variables 
(coefficients). For example, the radius of a circular waveguide can't be less than or 
equal to zero. Similarly there are limits to which n can be realistically varied for a 
dielectric. Routines for constrained optimization are available in standard C and 
Fortran libraries, and also in software packages like MATLAB available from The 
10 Math works, Inc. of Natick, Massachusetts. Most of these routines work best when 
the number of optimization variables is small, say 2. Use the hierarchical, sectional 
optimization procedure explained in the previous section to optimize P 0 as a 

function of the properties of the obstructing structure. As optimization is continued 
the number of variables may be increased to more finely define the obstruction, thus 

15 achieving more flexibility in optimization and reaching a desired value of Pq as a 
percentage of the total input power. 

If there is more than one parameter that specifies the output field, e.g., powers 
and phases of more than one mode, then the optimization becomes even more 
difficult. For example, consider the case where the output electromagnetic field has 

20 two modes with equal power in each. Then we can maximize the following two 
parameters: 

1. Po. the total additive power in the two modes, i.e., 

Po = Pi + where Pi and P 2 are the powers in the two modes 

respectively. 

25 2. R = Psmaii^Piarge* where Pgmau is the smaller of the two values Pj and P 2 

and Piarge iS * e larg er of the two values. 

Thus when carrying out the design of a specific mode 
control/conversion element, the additional items (for example) would need to be 
specified: 

30 

Output Field: 50% TE n and 50% TM n 

Optimize: Pq = Total output power in TEu(Pi) and TMn(P 2 ) 

modes. P sma ii is the smaller of the two mode 
powers Pi and P 2 and P^ge is the larger of the two 
35 mode powers. 

-n- 
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Optimize: R = Psmall/Plarge 

9. OBTAIN MULTIPLE DESIGNS (OPTIONAL) 
This step is optional but preferred. Repeat steps 5 through 8 to come up with 

a number of different designs based upon different initial starting points for the 
structures (step 6). Select one of these designs to build the mode control/conversion 
device based on simplicity of structure, bandwidth of operation, size, ease of 
fabrication, and the ability of the design to tolerate fabrication errors (maximum 
possible in the fabrication technique) without drastically affecting the value of Pq. 

10. FABRICATE THE DEVICE 

Fabricate the mode control device using an appropriate fabrication technique. 

DESIGN EXAMPLES 

15 What follows are examples which demonstrate the techniques used to design 

different types of mode control/conversion devices (in different waveguide 
geometries) using the methods of the present invention. The following descriptions 
are based upon the same basic steps described in the previous section. 

20 EXAMPLE 1: TE mn to TE pq Mode Converter in a Rectangular Waveguide 

Using Only Material Variation 

1. Specify the application of the device: 

(a) Frequency: 4 GHz 

25 Input Waveguide: Dielectric filled (n ) metallic wall 

rectangular 

waveguide with dimensions "a" in the x- 
direction and "b" in the y-direction. 
30 Output Waveguide: Same waveguide as at the input. 

(b) Input EM Field: A known combination of TE mn modes. 

(c) Output EM Field: Only TE^ mode. 

(d) Output Parameter Power in TEpq mode. 



-I*- 
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2. Pick a method of variation: 

Here it is desired to develop the mode converter by putting a dielectric 
material inside the waveguide 20 at box 22 as shown in FIG. 4. Thus, the method of 
variation in the material properties of the obstruction is selected and the waveguide 
dimensions remain constant. As the obstructing device 22 is a dielectric material, its 
magnetic permeability, p., is also assumed to be constant 

3. Choose the material: 

Choose a suitable material for the structure. The material should be such that 
its n can be varied significantly and various shapes can be fabricated easily using 
this material. A wide variety of such materials are known to those of ordinary skill in 
the art. Examples include: teflon and ceramic. 

4. Decide the direction of perturbation: 

Because it is desired to convert from TE mn modes to a TE W mode, i.e., both 
suffixes need to be converted, n will thus be allowed to vary in all the three 
dimensions x, y, and z. 

5. Pick a suitable basis function: 

The variation in n is represented by using sinusoidal basis functions such 

that: 

An (x,y,z) = x 0 sin(x) + Xxsin(3x) + ... + yosin(y) + yisin(3y) + ... + 
ZQsin(z) + zjsin(3z) + ... 

The series is truncated to two terms in each dimension, yielding: 

An (x,y,z) = x osin(x) + xxsin(3x) + y 0 sin(y) + yjsin(3y) 
+ ZQsin(z) + Z!sin(3z) 

6. Choose an initial structure approximation: 

Assign suitable values to XQ,x l? yo,yi,zo,zi and also chose lengths Lj, L2, and 

L 3 - 
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7. Formulate the forward solution: 

The forward solution is developed by writing a FEM (Finite Element Code) or 
by using HFSS software to find power P 0 in the mode TEpq. Then 

Po = P 0(*0' X I> y0> Vh Zfr z l> L l> ^ L 3> 

5 

8. Perform a global optimization: 

BCPOL, an IMSL (International Mathematical and Statistical Library) 
constrained optimization routine is preferably used for optimizing P 0 . Variables Xq, 

x i* yo» Yi» z 0» z i» L 2 and L3 are optimized using the hierarchical, sectional 
10 optimization procedure explained in the previous section. The optimized variables 
then define the structure of the dielectric obstruction for mode conversion. BCPOL 
is available from IMSL of 2500 City West Blvd., Houston, Texas. 

9. Obtain multiple designs (optional): 

15 Optionally repeat steps 5 through 8 to generate a variety of designs based 

upon different sets of initial values for the variables: xo,x 1 ,yo»yi,Z0,z 1 ,Li,L 2 and 

L 3 . Select the best design based on performance and ease of fabrication. 

10. Fabrication: 

20 Fabricate the final design using an appropriate technique. 

EXAMPLE 2: Free Space Optical Wavelength Filter : Bragg Reflector 

This example provides a methodology for preparing mode control/conversion 
25 devices for use in optical waveguide systems, e.g., optical fiber. In order to be able 
to more fully utilize the available optical fiber spectrum for communication of data, 
wavelength division multiplexing of the signal can be used. The limitation of the 
system becomes how close the lines can be spaced. The LASER line needs to be of 
very narrow bandwidth and stable in order for close spacing. In addition, narrow 
30 spectral linewidth (and hence long coherence length) light is important for 
metrology. 

Conventional periodic Bragg reflectors cannot produce as narrow a passband 
as aperiodic structures of the same total thickness. Rather than relying on end 
mirrors, the spectral and mode control structure could be distributed throughout a 
35 LASER diode cavity. Using lithography or special semiconductor material growth 
procedures known to those of ordinary skill in the art, it is possible to have materials 
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with different refractive index properties distributed throughout the LASER cavity. 
By controlling the size and location of these materials, it is possible to also control 
the optical signal properties of the cavity in a superior fashion to the standard cavity 
approach. The location and size of the "inhomogeneities" that are distributed 
5 throughout the cavity would be optimized to achieve the desired result. 

In a wavelength division multiplexing system, filters for multiple wavelengths 
are needed where several wavelengths are passed and the remainder of the spectrum 
is rejected. Standard filter theory can be used to design a bandpass filter for a single 
wavelength. If there are, for example, three wavelengths, then superimposing or 

10 cascading designs for each wavelength will not work very well because a filter for 
one wavelength will not pass the other two effectively. By using optimization 
techniques with aperiodic structures, better filters for the three wavelengths can be 
designed. These filters are implemented by changing the dielectric properties or 
profile of an optical waveguide (e.g., planar or optical fiber). 

15 Here, the goal is to design a narrow band filter for optical wavelengths to be 

used in narrow band distributer Bragg reflector ("DBR") or distributed feedback 
("DFB") LASERs. Currendy periodic dielectric variations, i.e., gratings or Bragg 
reflectors, are used in DBR LASERs to make them narrow band. In accordance with 
the technique of the present invention, these are replaced with irregular structures 

20 the shapes of which are not pre-defined. These structures generally achieve 
narrower bandwidth of operation than can be achieved by Bragg reflectors. Thus in 
this way the present invention can be used to design narrow wavelength filters for 
optical applications. 

25 1. Specify the application of the device: 



(a) 



Frequency: 



Optical LASER Frequency (Wavelength = 
1480nm ) 



(b) 



Input EM Field: 



Normally incident single plane wave in free 
space (say output of a LASER). 



30 



(c) 



Output EM Field: 



Normally exiting single plane wave in free 



(d) 



Output Parameter: 



space. 

Power in the output wave. 



- 2J 
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2. Pick a method of variation: 

For obvious reasons choose method 2b: Variation of material properties of 
the obstruction. 

3 . Choose the material: 

The scattering structure will consist of a thick dielectric slab 24 schematically 
shown in FIG. 5 whose dielectric constant varies in the z direction only (the 
direction of the incoming wave). It is preferable to use a material that is easy to use 
and can have maximum variation of n . Al x Gai_ x As is suited for such an 
application. Other suitable materials will also be readily apparent to those of 
ordinary skill in the art. By changing "x" in the chemical formula, a large variation in 
n can be achieved (8.8 to 12.8). Using an MBE (Molecular Beam Epitaxy) 
machine it is possible to grow layers of Al x Ga loc As with varying thicknesses. A 
structure as shown in FIG. 5 made up of slabs of Al x Gai. x As which have essentially 

infinite transverse dimensions in comparison to the spot radius of the input LASER 
beam is used to implement the filter. 

4. Decide the direction of perturbation: 

The dielectric variation will only be in the z direction. 

5. Pick a suitable basis function: 

For the type of structure proposed in step 3, an ideally suited basis function is 
the step function^ Thus: 

An (z) = eoU(z) + e 1 U(z-t 1 ) + e2U(z-t 2 ) + ... (U is the step function) 

where each of the variables eo» e2, ... can vary from 7.8 to 11.8 and An is the 
difference between the n of the material and the n of free space. The only limits 
on t lf t 2 , t 3 , ... are as follows: 

ti+i-ti^W, for alii 

ti+rti^W for alii 

t^: Minimum sheet thickness that can be generated by the MBE process. It is 

equal to approximately one atomic layer, 
tj^: Practically feasible maximum thickness of a slab. It is generally equal to a 

wavelength at the operating frequency of the LASER. 
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Truncating the series to a finite solution: 
p o = PoC^ e l' e 2—- » ^ l l> l 2* > ^ L ) 

5 

6. Choose as initial structure approximation: 

As an initial "guess" one may alternately pick values of ej's to be 7.8 and 11.8, 
select t i+1 -ti = 0.5 wavelength for all i's, and pick an initial length, L, of 5 
wavelengths. 

10 The dielectric constants (10 in number) and widths of the slabs (10 in number) 

are then the finite number of variables which uniquely identify the obstructing 
structure 24. 

7. Formulate the forward solution: 

15 The forward problem is easy to solve. There is a standard analytic expression 

available to solve for the transmitted power through a series of dielectric slabs. 

The transmitted power for a stack of three dielectrics is calculated as follows. 
See FIG. 5 A. Let n lf n 2* n 3 be the characteristic impedances in the three 
mediums 25a, 25b and 25c. Then the impedance seen at the input to medium 25b is 

20 given as: 

- n 2 (( n 3 cos ^2* + J n 2 s " 1 ^V( n 2 cos ^2* + J n 3 s * n ^2*)) 

where k2 is the wave number in medium 25b. Then the reflection coefficient is 

25 

given as: 

p = ((Z L -n !)/(Z L + n 0) 
^transmitted = ^incident 0~'P' 2 )- 

30 The output power, P 0 , at wavelength 1480nm is then a function of 

Pq^' e l» e 2»— ^txu *i» h* *nt L). 



-z3 - 
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8. Perform a global optimization: 

The FORTRAN IMSL routine BCPOL is preferably used to globally optimize 
P 0 . For this simple problem BCPOL could easily optimize up to 8 variables 
simultaneously. The constrained optimization limits given in step 5 for various 
5 variables are given as input to the BCPOL program. The sectional, hierarchical 
optimization technique is then used to determine the optimized structure. As 
optimization progresses, the number of optimization variables is increased for more 
flexibility by adding more slabs at the beginning or end of the structure, thus 
increasing L. The process is stopped when almost 100% power is transmitted 
10 through the structure at the nominal wavelength. 

9. Obtain multiple designs (optional): 

Using different initial guesses and design decisions (as in step 6) a variety of 
structures are designed. The 3 dB bandwidth for each structure is obtained by 
15 solving the forward electromagnetic problem as a function of frequency, i.e., by 
evaluating the transmitted power at various frequencies and plotting power versus 
frequency to determine the 3db points. The device with the least bandwidth is 
chosen for fabrication. 

20 10. Fabricate the device: 

The dielectric slab structure is fabricated by growing layers of 
Al x Gaj_ x As by MBE (Molecular Beam Epitaxy) or any other suitable technique. 

EXAMPLE 3: Circular Waveguide TE n to TM n Microwave Mode Converter 
25 from Smaller to a Larger Waveguide 

This structure can be called a mode converter, a mode launcher, or a 
waveguide transition. It is a mode converter because it converts a TEn niode to a 

TMn mode. It is a mode launcher because it demonstrates that any type of mode 
30 can be launched into an overmoded large waveguide which is fed by a small 
waveguide having a single known mode in it. 

1. Specify the application of the device: 



(a) 



Frequency : 
Input Waveguide: 



9.94 GHz (n =3 cm). 

Air-filled metallic wall circular waveguide 
with radius 0.3895". 
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Output Waveguide: Air-filled metallic wall circular waveguide with 



radius 0.812". 



5 



(b) 



(c) 



(d) 



Input EM Field: Only composed of TEj j mode. 

Output EM Field: Only composed of TMn mode. 

Output Parameter Power in TMj j mode. 



2. Pick a method of variation: 

Method 2a is selected (variation of waveguide shape) because this is easy to 
demonstrate and also results in a new way of fabricating microwave waveguide 
10 mode converters, as will be discussed later. 

3. Choose the Material: 



15 4. Decide the direction of perturbation: 

Because both TEn and TM n modes are circularly symmetric, a change in the 

radial dimension of the waveguide as a function of z would be suitable for such a 
conversion. 

20 5. Pick a suitable basis function: 

In order to demonstrate a very inexpensive and modular fabrication 
technique for microwave mode converters select step functions to represent the 
delta variation in radius as a function of z. Then: 

25 Ar(z) = r 0 U(z) + rjLKz-Z!) + r 2 U(z-z 2 ) + ... . (U is the step function) 

The following restrictions are used on the range of possible values for 
ro,r lt r 2 , ... and zj's: 

ri's should be selected so that the radius of the mode converter section 
30 at no place is zero or less and is not greater than r max (1.75"). 

z i+i " z i = As where the value of As, the step size, will be a variable in 
the optimization procedure, i.e., the power Pq in the TM n mode will depend upon 
the following values: 



To keep it simple, choose an air-filled mode converter section. 



35 



Po = p 0 < r 0> r i> r 2 r n. A*, L) 
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where the series has been truncated to n+1 values. 

Note that in this formulation positive values of r^s are permitted, i.e., the 
structure can protrude outside the walls of the waveguide (see, e.g., FIG. 3). It has 
5 been found through experimentation that protrusions outside the waveguide results 
in quicker convergence of the optimization solution and also results in shorter mode 
converters. What should be the size of the largest protrusion is a design decision. 
As the largest protrusion radius, r max , is increased, the computation time for the 

forward solution increases exponentially. In some cases the computation time could 
10 be prohibitive to achieve optimization convergence. Also it is desirable to limit the 
radial size of the mode converter to keep the structure compact. 

6. Choose an initial structure approximation: 

Initially only four terms from the series of step functions are included and use 
15 the following values: 

ri = 0.2355", r 2 = 0.25", r 3 = 0.375", r 4 = -0.25", 

As = 0.2505" and L = 1.002" 

It is desired to fabricate sections of this mode converter out of commonly 
20 available 0.2505" thick aluminum sheet. Therefore, the As are held constant and are 
not included as a variable in the optimization procedure. This is not required as the 
As could be varied as well, if desired. 

7. Formulate the forward solution: 

25 The mode converter structure defined above consists of equally spaced 

waveguide discontinuities. Given the information in steps 1 through 6 the power in 
the modes generated by the structure can be calculated using a variety of methods. 
Two of these methods are Finite Element Method (FEM) and the Mode Matching 
Method. Hewlett Packard's commercial software HFSS also has the capability to 

30 solve for waveguide discontinuities. The power P 0 can therefore be written as a 

function of the various variables, i.e., 

Po = p o( r O> r!,r 2 ,r 3 ,L) 

8. Perform a global optimization: 

35 BCPOL is again preferably used for constrained optimization of the variables. 

Constraints as given in step 5 are used as inputs to the BCPOL algorithm. It was 
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found that the BCPOL solution, in this case, converged faster if the number of 
variables were restricted to 2. Thus, sectional optimization was employed. To 
reduce the total time in the calculation of the forward solution, the scattering effect 
of the whole structure is also calculated in sections. While two of the variables are 
5 being optimized, the solution for the other variables is only calculated once and 
stored in memory. This technique saves a great deal of time, especially when there 
are a large number of variables. 

In the hierarchical optimization procedure, the number of variables can be 
increased in two ways; by reducing the value of As, or by increasing the overall 
10 length of the structure. The technique of reduction in As can be used until As^ is 

reached, which is the minimum acceptable thickness of each disc in the structure. 
As m j n is chosen as a fabrication parameter. In this case As = As m i n = 0.2505'*. Thus, 

the only option that there is to vary is to increase L (i.e., the number of discs). The 
value of L is increased gradually during the hierarchical procedure, which increases ^ 
15 the domain of possible structures. The procedure is stopped when power in the 
TM n mode is close to 100%. 

9. Obtain multiple designs (optional): 

Using different initial guesses and design decisions (as in step 6; when and 
20 where to add more steps etc.) various mode converters were designed. The' 
bandwidth of operation was calculated for each design and also estimates were ' 
generated as to how much the operation of each was affected by incorporating 
realistic fabrication errors. Based on these factors the best design was selected for 
construction. The total length of the finally designed mode converter is 2.505" = 
25 10*0.2505" and its conversion efficiency is 99.5%. A cross-sectional view of the 
TEji to TMn circular waveguide mode converter is shown in FIG. 6. 

10. Fabricate the device: 

The waveguide mode converter may preferably be fabricated as described 

30 below. 

For microwave mode converters one of the biggest advantages of the present 
invention is that the designed mode converters are very compact and short in size 
compared to earlier designs based on periodic wall perturbations (grating structures). 
There are two design parameters that affect the overall length of the designed mode 
35 converter, As min and r max . By reducing As min and by increasing r max the overall 

length of the mode converter can be reduced. There are, however, practical limits on 

-27- 
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the minimum value of As min and the maximum value of r max . The value of As m j n is 
restricted by the fabrication procedure, material used, and fabrication errors. 
Whereas the maximum value of r max is restricted by the available computation time 
and physical restrictions on the radial size of the structure depending upon its usage. 
5 Instead of using earlier techniques of pre-forming or molding for the 

fabrication of microwave mode converters, as with the periodic structures discussed 
above, a new modular approach to the construction of waveguide mode converters 
may preferably be used. In this approach the structure of the mode converter is 
discretized along the z dimension as shown in FIG. 6. Each discretized section is 

10 then like a disc which can be machined individually. All the discretized sections can 
then be put together inside a housing to assemble the complete structure. It was to 
demonstrate this idea that step functions were used in the design procedure for the 
TEn toTMu mode converter. Continuous basis functions could have also been 
used in the design (e.g., sine functions). In that case the continuous structure of the 

15 mode converter would have to be discretized to implement the modular fabrication 
method. In such a discretization the difference in the radius of two consecutive 
discs and the width of every disc will have to be kept less than about one-twentieth 
of a wavelength. This is required to maintain an accurate representation of the 
continuous surface by using a very small step size compared to a wavelength. This 

20 will increase the number of discs to be machined. To keep the demonstration of the 
idea simple, discontinuous basis functions have been used which reduce the number 
of required discs in the structure. 

The assembly of the mode converter 26 is shown in an exploded view in FIG. 
7. The thickness of each disc 28a, 28b, 28c, 28d, 28e, 28f, 28g, 28h, 28i, 28j in the 

25 preferred mode converter of FIG. 6 is 0.2505" and there are a total of 10 discs. An 
aluminum sheet with thickness slightly greater than 0.25" is machined on both sides 
to smooth its surfaces and to have a uniform thickness of 0.2505". A precision 
grinding machine is used for this purpose. Ten squares with side dimension slightly 
greater than 4.44" are cut out of the aluminum sheet. The square sections are then 

30 machined to make them into circular discs, each of which has a outer diameter of 
4.44". This is the ID (inner diameter) of a 4.5" OD (outer diameter) aluminum tube 32 
that will house these discs for concentricity. Appropriate sized holes 30a, 30b, 30c, 
30d, 30e, 30f, 30g, 30h, 30i, 30j are drilled in the center of each square (the radii of 
these holes are the same as the optimized radii of discs in the mode converter 

35 design). 
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The ten discs are press fit into aluminum tube 32. The length of aluminum 
tube 32 is about 1.5" larger than the total length of the mode converter 26 (2.505"). 
This space is to slide in the input and output waveguide support units 34, 36, 
respectively. These support units are machined out of thick aluminum pieces. The 
input unit 34 has a center hole 38 large enough for the input waveguide 40 while 
the output unit 36 has a hole 42 to slide in the output waveguide 44. The input and 
output adapters have vertical set screws 46, 48 to hold the two waveguides 40, 44 
in place. The input and output waveguides 40, 44 are preferably pieces of 
commercially available copper and brass tubing with average OD 0.875" and 1.75", 
respectively, and with wall thicknesses of 0.032" and 0.063", respectively. The 
whole structure of the mode converter is assembled together by a set of four 
through holes, nuts and bolts as shown in FIG. 7. The fully assembled mode 
converter 26 is shown in FIG. 8. 

The radii of the input and output waveguides 40, 44, respectively, are 
0.3895" and 0.8120". The radii of the holes in discs 28a - 28j for this particular 
application are, respectively, 0.3543", 1.2874", 1.5630", 1.5748", 0.8268", 0.9055", 
0.8268", 0.9055", 0.9449" and 0.7756". 

All components used in this design were commercially available and no 
special machines were used in the procedure. Due to this reason the cost of 
fabrication was very low, unlike more conventional mode converters. 

In order to test the performance of the mode converter a TE n mode was first 
excited in the input waveguide 40. A microwave source having a coaxial cable 
output was used for this purpose. An extended center-pin SMA launcher shown in 
FIGS. 9-10 was mounted on the side wall of the input waveguide of the mode 
converter. The radius of the input waveguide is such that at 9.94 GHz it only 
supports the dominant TE n mode. The SMA launcher, due to its impedance 

mismatch with the waveguide, not only launches the TE X1 mode but also a number 
of higher order evanescent modes into the waveguide. The length of the input 
waveguide from the SMA launcher to the input plane of the mode converter is 
selected to be large enough so that any evanescent modes have negligible amplitude 
at the mode converter plane. Thus is ensured a pure TE U mode at the input of the 
mode converter. The output of the mode converter was tested using a far field 
measurement technique and was found to give a 98.1% conversion efficiency to 
TM n mode, which is very close to the calculated value of 99.5%. 
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DISCUSSION OF INITIAL STRUCTURE APPROXIMATION 



It is logical to ask how arbitrary and wrong can the initial guess be and still 
have the process work? A solution will be found regardless of the initial guess for 
the iterative optimization solution. How good the solution is depends on the type of 
5 optimization and the closeness of the initial guess to the final, good solution. With 
the optimization solution in a mode converter application, the goal is to maximize the 
power in the desired output mode, thereby minimizing the power in all other 
reflected and transmitted modes. Another way to view this is that the wasted power 
is being minimized. The optimization surface, which is a function of the variables 

10 used (such as waveguide dimension), is very rough, meaning that there are many 
local minima, i.e., valleys whose minima are larger than others (which may be close to 
the desired zero wasted power and hence global solution). The goal of an 
optimization routine would be to find the global minimum by jumping out of the 
local valleys in search of a global solution, meaning the best given the parameter 

15 search space. If the optimization algorithm is very good, any initial guess can be 
used. Many approaches are available to accomplish this task in addition to those 
described herein. 

FIG. 1 1 shows an example of a simple scattering body 50 within a waveguide 
52 of diameter 15cm at a center frequency of 10 GHz. Scattering body 50 is formed 
20 of two 0.5 n -thick sections, one of height h] and the other of height h 2 - The z 

dimension of each is the same. In FIG. 12 is shown the optimization surface for TEj 

to TE2 conversion for the situation depicted in FIG. 1 1 . This is a very simple case. 

As can be seen easily in FIG. 12, the wasted power is depicted as increasing along 
the vertical axis in the upward direction. The goal of optimization of Pq, for example, 

25 is to obtain a solution with a low value of wasted power (optimized Pq). Thus the 
"lower" points on the optimization surface which correspond to the indicated step 
heights for h] and h 2 in units of 0.1 n are "better" places to be than the local 
maxima of wasted power. While a "best" solution may be impossible or difficult to 
achieve under many circumstances, many good solutions are available. Thus one 

30 reason for pursuing several different starting points is to ensure that one of the 
"good" solutions will be obtained, rather than a mere local minimum of the surface. 

FIG. 13 depicts an optimization surface for TEj to TE$ mode conversion 

under the circumstances of FIG. 11. As can be seen, the surface here is much more 
complex than that of FIG. 12 with many more local minima demonstrating the 
35 desirability of starting from a variety of initial points to guarantee a "good" solution. 



3o- 



WO 98/32187 



PCT/US98/00523 



ALTERNATIVE EMBODIMENTS 
A number of different solutions to mode conversion and control have been 
obtained to date. Several of these are shown in the figures. 

FIG. 14 is a diagram showing the construction of a parallel plate waveguide 
TEj to TE 2 converter in accordance with a presently preferred embodiment of the 
present invention. 

FIG. 15 is a series of four diagrams showing the sequence of optimization of 
the structure for a parallel plate waveguide TEi to TE 2 mode converter structure in 
accordance with a presently preferred embodiment of the present invention. 

FIG. 16 is a diagram comparing three distinct designs for a parallel plate 
waveguide TE j to TE 2 mode converter structure in accordance with a presently 
preferred embodiment of the present invention. 

FIG. 17 is a chart comparing performance parameters of the three distinct 
designs shown in FIG. 16. 

FIG. 18 is a diagram showing the structure of a parallel plate waveguide TEx 
to TE 2 mode converter structure in accordance with a presently preferred 
embodiment of the present invention. 

FIG. 19 is a diagram showing the structure of a circular waveguide TEq 2 to 
TEqi mode converter structure in accordance with a presently preferred embodiment 
of the present invention. 

FIG. 20 is a plot of efficiency versus frequency for the circular waveguide 
TE 0 2 to TE 0 i mode converter structure of FIG. 19. 

FIG. 21 A is a plot of width profile versus distance along the z direction for a 
prior art circular waveguide TEq 2 to TE 0 i mode converter in accordance with M. /. 
Buckley, et aL, "A Single Period TE 02 to TE 0J mode converter in a highly 
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overmoded circular waveguide," LE.E.E. Trans. MTT, Vol. 39, No. 8 (Aug. 1991), 
pp. 1301 - 1306. 

FIG. 2 IB is a plot of efficiency versus frequency (GHz) for a prior art circular 
5 waveguide TE02 to TEqi mode converter in accordance with M. 7. Buckley, et aL, "A 
Single Period TE 02 to TE 0l mode converter in a highly overmoded circular 
waveguide, " LE.E.E. Trans. MTT, Vol. 39, No. 8 (Aug. 1991), pp. 1301 - 1306. 

FIG. 22 is a chart comparing the length and conversion efficiency of the 
10 "Previous" design of FIGS. 21 A and 21B with the "Present" design of FIGS. 19 and 
20. 

FIG. 23 is a diagram of a grating design for a circular waveguide TE05 to TEqi 
mode converter for operation at 140 GHz in accordance with a presently preferred 
15 embodiment of the present invention. 

FIG. 24 is a plot of efficiency versus frequency (GHz) for a prior art circular 
waveguide TE 06 to TE^ mode converter for operation at 140 GHz in accordance 

with K. Kumric et aL, "Optimization of Mode Converters for Generating the 
20 Fundamental TE 0I Mode from TE 06 Gyrotron Output at 140 GHz, " International 

Journal of Electronics, Vol. 64 (Jan. 1988), pp. 77 - 94. 

FIG. 25 is a chart comparing the length and conversion efficiency of the 
"Previous" design of FIG. 24 with the "Present" design of FIG. 23. 

25 

FIG. 26 is a pair of diagrams comparing a pair of designs for a circular 
waveguide TE n to TMn mode converter structure each designed in accordance 
with a presently preferred embodiment of the present invention. 

30 FIG. 27 is a chart comparing the length and conversion efficiency and other 

performance parameters of the "Design 1" design of FIG. 26 with the "Design 2" 
design of FIG. 26. 

Although illustrative presently preferred embodiments and applications of this 
35 invention are shown and described herein, many variations and modifications are 
possible which remain within the concept, scope, and spirit of the invention, and 
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these variations would become clear to those of skill in the art after perusal of this 
application. The invention, therefore, is not to be limited except in the spirit of the 
appended claims. 
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CLAIMS 

What is claimed is: 

1. A method for creating a design for a waveguide mode control device 
5 comprising the steps of: 

determining the application of the mode control device including its 
frequency of operation, type, structure, size of the input and output waveguides, 
mode composition of the incident electromagnetic field and the required mode 
composition of the output electromagnetic field; 
10 determining a method of variation of an obstruction to be placed 

within the mode control device; 

choosing a material composition of said obstruction; 
determining the directions in which said obstruction may vary; 
determining a suitable basis function to represent the variations in the 
15 obstruction and form an expression utilizing combinations of variables and said basis 
function to represent variation of said obstruction in at least one dimension; 

selecting an initial structure approximation of the obstruction and 
thereby selecting actual initial values for said variables and a length for the mode 
control device; 

20 formulate a solution for output power in at least one selected output 

mode as a function of said variables and said length; 

perform a global optimization to determine a set of values of said 
variables and length that tend to maximize output power in said at least one selected 
output mode. 



25 



30 



2. A method according to claim 1 , further comprising the step of: 
repeatedly carrying out steps f, g and h with different sets of said initial 

values to determine different sets of values of said variables and length that tend to 
maximize output power in said at least one selected output mode. 

3. A method according to claim 1 , further comprising the step of: 
fabricating the mode control device in accordance with said set of 

values of said variables and length that tend to maximize output power in said at 
least one selected output mode. 



35 



34- 



WO 98/32187 PCT/US98/00523 

4. A method according to claim 2, further comprising the step of: 
fabricating the mode control device in accordance with one of said 

different sets of values of said variables and length that is nearer a desired solution 
than others of said different sets of values. 

5 

5. A mode control device for use with a first optical waveguide and a 
second optical waveguide, said mode control device comprising: 

a first planar layer of material; 

a second planar layer of material adjacent to said first planar layer of 

10 material, 

said first planar layer of material disposed so that electromagnetic 
energy from said first optical waveguide passes through said first planar layer of 
material and then through said second planar layer of material and eventually into 
the second optical waveguide, 
15 said first planar layer of material having a first index of refraction; 

said second planar layer of material having a second index of 

refraction. 

6. A mode control device according to claim 5, wherein: 

20 said first planar layer of material comprises Al x Ga!_ x As for a first value 

of x. 

7. A mode control device according to claim 6, wherein: 
said second planar layer of material comprises A^Ga^As 

25 for a second value of x, said second value of x being different from said first value of 
x. 

8. A method according to claim 1 , wherein: 
said basis function is a delta function. 

30 

9. A method according to claim 1 , wherein: 

said basis function is the Kronecker delta function. 

10. A method according to claim 1 , wherein: 
35 said basis function is a step function. 



WO 98/32187 



PCT/US98/00523 



11. A method according to claim 1 , wherein: 
said basis function is a Bessel function. 

12. A method according to claim 1, wherein: 

said basis function is a trigonometric function. 

13. A method according to claim 12, wherein: 
said trigonometric function is a sine function. 

14. A mode control device for use with a first circular microwave 
waveguide supporting a first set of modes and a second circular microwave 
waveguide supporting a second set of modes, said mode control device comprising: 

a first disk having a first aperture of a first radius; and 

a second disk adjacent said first disk having a second aperture of a 
second radius, said second radius being different from said first radius, 

said first disk disposed so that electromagnetic energy from said first 
circular microwave waveguide passes through said first aperture and then through 
said second aperture and eventually into the second circular microwave waveguide, 

microwave energy within said first set of modes incident upon said first 
aperture scattered into said second set of modes in said second circular waveguide 
and into said second aperture, 

microwave energy within said first set of modes incident upon said 
second aperture scattered into said second set of modes in said second circular 
waveguide. 
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